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We recently generalized in the analytical Bethe Ansatz for integrable open spin chains with 
quantum- algebra invariance which was developed in to the entire series of C/g(i?„ )-invariant 
spin chains in the fundamental representation. (The analytical Bethe Ansatz for closed spin chains 
with periodic boundary conditions was formulated by Reshetikhin We focused in |T[ on the 

(2) 

A>2^ series because we sought to identify the main difficulties in generalizing the analytical Bethe 
Ansatz procedure to any affine Lie algebra, and the series was particularly convenient since 

the Izergin-Korepin Q ^2 case was already understood Q. The main difficulties were computing 
the pseudovacuum eigenvalue of the transfer matrix, and formulating an appropriate Ansatz for 
general eigenvalues. 

In dH a "doubling postulate" (i.e, that the Bethe Ansatz equations are "doubled" with respect to 
those of the corresponding closed chain with periodic boundary conditions) was introduced. Using 
this "doubling postulate", we were able to easily formulate an appropriate Ansatz and obtain the 
transfer matrix eigenvalues. [| Very recently this procedure was used for the spin chain by 
Yung and Batchelor |Q. These authors have further generalized this method to certain open spin 
chains which are not quantum-algebra- invariant ^, |6|. 

The eigenvalues of the transfer matrix have been obtained (by means of the algebraic Bethe 
Ansatz) also for the An ^ open spin chains 0, |8|, and the B^^^ open spin chain 

The success of the analytical Bethe Ansatz procedure gives us confidence that the same pro- 
cedure should work for the remaining series of quantum- algebra- invariant open spin chains, with 
Hamiltonian H = ^Yl!j=\ I '^'^\u=q ■ Here K{u) is the R matrix associated with an affine 

Lie algebra g^^^ and with the fundamental representation of 50 1 where go C g^^^ is the maximal 
finite-dimensional subalgebra of g'^^^ ■ Unfortunately, the (diagonal) K matrix which is needed to 
construct the corresponding transfer matrix is known only for the following additional series 
of R matrices: J^2n-\^ Bn \ Cn \ Dn^ . We therefore restrict ourselves to these cases. 

Specifically, in this paper we determine the eigenvalues of the transfer matrices for the following 
four infinite series of quantum-algebra invariant open spin chains: 

• C/q(C„)-invariant A2n-i ^V^^ chains 

• C/q(i?„)-invariant Bn^ (n > 1) spin chains 

• C/q(C„)-invariant Cn^ spin chains 

• C/q(D„)-invariant Dn^ spin chains 

^This postulate provides a short-cut for obtaining the transfer matrix eigenvalues. In principle, this postulate 
can be avoided, and the transfer matrix eigenvalues can be obtained by carefully implementing the constraints of 
analyticity, crossing, fusion, asymptotic behavior, and periodicity. In practice, however, this can be quite tedious. 
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We emphasize that within the class of integrable open spin chains, those spin chains which 
are quantum-algebra invariant are in certain respects the simplest ones. Indeed, for the quantum- 
algebra-invariant spin chains, the Bethe Ansatz states are highest-weight states of the quantum 
algebra |11] - |14]. Moreover, for those open spin chains which are not quantum-algebra invariant, 
the Bethe Ansatz equations have additional factors depending on some boundary parameters. 



We therefore work with solutions R{u) [T^, 16 1 of the Yang-Baxter equations which are associated 
with the affine Lie algebras g^'^^ = (^2n-i i^n^ ,Cn^ ,Dn^) and with the fundamental represen- 
tation of the Lie algebras qq = {Cn , Bn , Cn , Dn) , respectively. The corresponding matrices R{u) 
commute with the generators of the quantum algebras Ug{go)- In the appendix we collect the 
necessary information about these solutions. We follow the notations of 
Our goal is to determine the eigenvalues of the transfer matrix 0, IC] 

t{u)=tVaMaTa{u)fa{u), (1) 

where 

Ta{u) = RaN{u) RaN~l{u) ■ ■ ■ Ral{u) , 

fa{u) = Rla{u) ■ ■ ■ RN~la{u) RNa{u) , (2) 

with the subscript a denoting the auxiliary space while the subscripts 1, • • • , refer to quantum 
spaces. The matrix M is given in the appendix. It is related to the crossing matrix V, M = V^ V, 
where 



Ri2{u) = Vi Ri2{-u - pY^ Vi (3) 



with 



.(2) 

-2Kr] for B^r!\ C'n> , D)t 



-i-K-lKT] for A\^_^ 



and K = (2n, 2n - 1, 2n + 2, 2n - 2) for B''n \ C^}\ D^n^), respectively. The transfer 

matrix commutes with Uq{go) | p!7| , |l2[. We consider simultaneous eigenstates of the transfer matrix 
t{u) and the n Cartan generators {Hi , • • • , Hn} of Uq{go). We call the corresponding eigenvalues 
^(mi ,m„) ^^-j ... ^ ^^1^ respectively. The eigenvalues of the Cartan generators are related 

to the integers mi , • • • , by^ 

{A/} = {N - mi ,mi - 1712 ■ ■ "in-i - "in} for b!^^ , 



^We correct a typo in M for the case C, 



(1) 
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{A;} = {iV - mi ,mi - 771-2 ,"i„_2 - m„_i - m„ ,m„_i - m„ } for D'^\ 
{Xl} = {N - mi ,mi - m2 ,■ ■ ■ , "^r^-l - 2m„} for , C^^ . 

We choose A^^''^ '"' '^"\u) to correspond to a highest weight vector for the corresponding algebra 
U,{go). 

To accomphsh the analytical Bethe Ansatz program we must specify the following additional 
information: 

1. Crossing relation M 



A(™i.-,mn)(^) = A(™1'-'™")(-u-/9), (5) 



with p given by Eq. 
2. Fusion formula H 



-C(n + /5)2^ 5(2tx + p) 5(-2tx-3p)}. (g) 



where 



C(n) = (7) 



= J 2sinh(| + 2,)cosh(H + .,) for 

^ 2sinh(f-2?7)sinh(^ + Kr/) for 5^,(7^^^^^ 



and 



COsh(| — KT]) for ^2n-l 



/3(n) = sinh(u - 2K7?) (9) 



sinh(| -K??) for B^n \ C^n \ D^n^ 



The fusion formula will be used in order to check the correctness of the pseudovacuum eigen- 
value. 

3. The transfer matrix is periodic t{u) = t{u + 27rz), and its eigenvalues A^*"! (ti) are 
analytic functions of u. 



3 



To obtain the pseudovacuum eigenvalue of the transfer matrix we compute its expectation value 
in the pseudovacuum state for = 2, 3 and we obtain 



C C ^mm \^ " /v'' ^mm) 

+ 62^{po(l + — ) + e^^^'^'^i - + — ]} 



rP fP- p2 j2 

+ a^^ eAWr?n , + ""^i , I /-^qn 

"mm e 2 _h2 ^ „2 - ^ (a? - r^Mo^ - b"^) ' ^ ' 

"'mm " "'mm ^ \"'mm ^ jV^mm " I 

where pq = I]^2^ wit)! Ma being the matrix elements of M given in Eq. (|^); d? = Y.'J^L^ a.^^; 
and a^p, b, c, and e are given in Eqs. ®, (IID in the appendix, m = 2n for ci^\ £>i^^ 

and m = 2n + 1 for Bn^ . 

We postulate that Eq. ( |lO|) is true for all A^. Using Mathematica we find the following expression 
for the pseudovacuum eigenvalue 

for 41,: 

A(0'-'0)(u) = c{uf^ sinh(n-2K7?) cosh(n - a;??) 

sinh(ii — 2rj) cosh(n — kt]) 



_ sinh(u) sinh(u — 2k?7) 



sinh(n — 2r/) sinh('u — 2(k — l)r/) 
, 2Ar sinh(n) cosh(n — {2k — uj)rj) 



+ sinh(.-2(.-l),)cosh(.-;,) 



for I^i^^ 



a(0,-,0)^a _ ^2jV sinh(u - 2/tr?) sinh(^ - a;r?) 

sinh(u — 2?7) sinh(u — Kr/j 
, , , sinh(u) sinhfu — 2Kri) 



sinh(n — 2ri) sinh(M — 2{k — 1)?]) 
2M sinh(u) sinh(ii — (2k — a;)r/) 



+ sinh(n-2(.-l),)sinh(n-;,) ^^^^ 



where 



2sinh(K??) cosh((K-2)?;) r r^(l) 
sinh(2»7) ' ^ 

Po = { (13) 

ln—2)ri) cc 
sinh(2r7) 



2smh((2«-2)>;)cosh(2nr?) r .(2) ^(1) 

lOl ^2n-l''-^" 
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_i K-2 for 

""I K + 2 for ^ ^ 

and c(u), b{u), amm{u) are defined in the appendix. This eigenvalue is consistent with the fusion 
equation (^), i.e. the expression for A^'^ '^^(n) does not have poles for values of u for which a{u) 
or /3(m) are zero. 

For the dressing of these eigenvalues we make the following Ansatz 
for 41,: 

sinh(ii — 2r]) cosh(n — Krf) 

+ C^"''\'-d)a„„n{"d?^ sinh(n) cosh(n - (2k - uj)ri) 

™™' sinh(ii — 2(k — l)r/) cosh(n — K'q) 

+ binf' I E + zKn) A^'"' | , (15) 

for Bi'\d'\D^P: 

sinh(ti — 2r]) sinh(n — kt]) 
n ^,^2W sinh(u) sinh(u - (2k - w)r?) 



sinh(ii — 2(k — I)//) sinh(ii — kt]) 

{n-l 

+z,(-) A^'"''"^'^^^(-)]}, (16) 



where 1=1 for and / = in the other cases. The functions A,B,C are the doubles of the 
corresponding expressions given by Reshetikhin Q (apart from slight changes in notation) and so 
are also invariant under u^P — > —uP 
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^(-i)(n) = Ji sinh(i(n - Q + t?) sinh(i (n + ^) + r?) ^ ^^^^ 
j=i sinh(i(n — u^p) — tf) sinh(i(ii + u'p) — rj) 

C(™i)(n) = A("i)(-u-/9) 



ijp''"'+^^(n) = n 

j=i smh{^ {u — u- ) — Irj) sinh(i(ii + uj- ) — /r/) 
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1 = 1, ••• ,n-2 for 4n-l,^^n^ 
/ = 1 , • • • , n - 1 for Si^^ 
/ = 1 , • • • , n - 3 for 
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for A! 



(2) . 
2n-l- 



for 



for a 



(1). 



(m„_i,m„) 



n-1 



for A 



(1). 



B, 



{mn-2,mn-l,mn) 

n-2 



= zi{-u-p), 



p), 



Z = 1, 



n 



sinh(i(n - nj""^^) - (n + l)r/) smh^u + nj""^^) - (n + I)?/) 



1 sinh(i(n — n^" ''"^) — (n — 1)77) sinh(^('u + n^" — (n — I)??) 



sinh(n 



2(n - 2)r?) smh(u + ttf ^ - 2(n - 2)r?) 



sinh(M 



,(n) 



2ra?7) sinh('u + u^"^ — 2n?7) 



n 



sinh(i(u - tij"^) - (n - 2)r/) sinh(^(u + uj"^) - (n - 2)77) 



sinh(^(M — Uy"^) — nri) sinh(^(u + u^J^^) — nrf) 



(n)^ 



sinh(i(ii 



in). 



(n + 1)77) smh(i(ti + uf^) — (n + I)??) 



sinh(i(« — 



n 



1)77) smh(i(t( + w^"'') — (n — 



n 



1 sinh(i(u - uj""^^) - (n + l)r/) sinh(^(ti + uj""^^) - (n + 1)77) 



1 smh.{\{u — Uj^ — [n — 1)77) sinh(^('u + — {n — 1)77) 



n 



sinh(i(u — u.!^' ) — (n — 3)7/) smh(i(u + Uj) — {n — 3)r]) 



1 sinh(i(u — Uj'^^) — {n+ 1)77) smh(i(7/ + u^"'') — (rz + I)??) 



sm]i(i(» - ''^) - ;?//) sm]i(i((/ + (/y' "^') 



I" I) 



(19) 



(20) 



(21) 



(22) 



J^i smh{^{u — uf ^^) — {n — 2)77) sinh(i(u + u^" ^'*) — (71 — 2)77) 
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"Jj^ sinh(i(u - uf^^O -{n- 3)7?) sinh(i(n + uf'^^) - (n - 3)r/) 
j=i sinh(^(u — -u^" ^^) — (n — 1)77) sinh(^(n + n^" "'^'*) — (?2 — l)r/) 
g smh(i(u - nj"^) - (n - 3)??) sinh(i(u + u^"^) - (n - 3)r?) 
j=i smh(^(u — ttj"^) — (n — l)??) sinh(^(u + Uj"^) — (n — I)//) 



^^smh(|(u — ) — (n — 3)7/) sinh(^(n + u^" ) — (n — 3)r/) 
=1 smh(i(u — u^" ^^) — (n — 1)7/) sinh(i(n + u^" ^^) — (?2 — l)r/) 
g sinh(i(n - uf'') -{n+ l)r/) smh(i(u + uf^) - (n + 1)??) 
j=i sinh(i(u — Wj"^) — (n — l)ry) smh(^(ti + ^"'*) — (n — l)r/) 



We notice that in contrast to the analytical Bethe Ansatz procedure for closed spin chains, Eqs. 
(|T5|) , (p!^) contain also the unknown functions zi{u) ,w{u). As mentioned in the introduction, these 
functions can be determined by the so called doubling postulate, i.e. we demand that the Bethe 
Ansatz equations obtained from the cancellation of poles in A, i3,C be "doubled" with respect to 
those in Reshetikhin's paper Q. The doubled Bethe Ansatz equations are 



sinh(-| r]) 

sinh(^^| — h r]) 



2N 



X 



7711 



= n 



sinh(i('U^^'' — Uj^^) — 2r]) sinh(i 



-rfcsinh(i(4^^ 



2r,) 



(1) 



m2 



n 



sinh(i( 



(1) (2) 
(2) 



+ 2r?) sinh(i(4^^ + n5-^^) + 2??) 
) + -q) sinh(i(4^^ + uf^) + r?) 



:i sinh(i(4^^ - u^p) - rj) sinh(i(4^'* + uf^) - i]) 



(24) 



sinh(l (4^ - uf^>) + r?) sinh(i(4^ + uf^>) + r?) 
j=i sinh(i(ii^^'* - Uj~^^) - v) sinh(i(4''' + - rj) 

^ sinh(i(4^ - 4^) - 2f?) smh(i(4^ + 4^) " 
j}k sinh(i(4^ - uf) + 27?) sinh(i(4'^ + 4^) + 2^?) 

sinh(i(4^ - 4""'^) + ^) sinh(i(4) + 4""'^) + v) 
j=i sinh(i(4^ - 4^^^) - V) sinh(i(4^ + 4^^^) ~ ^) ' 
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where 

Z = 1 , • • • , n - 1 for B'i^ 

l = l,--- ,n-2 for ^li^C-n'^ (26) 

Z = 1 , • • • , n - 3 for o'^n^ . 

Moreover, the Bethe Ansatz equations corresponding to values of I = 1 , 2 , • • • , n which are not 
included in Eq. ( p6[ ) are as follows: 

for B^n^: 



for 41i: 



"'-^ 5inh(i(4") - u'f-'^) + 7?) sinh(l(txi") + -^)) + r?) 
j=i sinh(i(n^"''* - u^""'^^) - ??) sinh(i(ti["'' + 4"^"'^^) - rj) 

sinh(l(^i") - nf )) - ^) sinh(i(4") + 4"^) - rj) 
,U sinh(l(4") - )) + ^) sinh(i(4") + 4")) + ry) 



"^-^ sinh(l(4"-^^ - 4"-^)) + r?) sinh(i(4"-^^ + 4"-^^) + r?) 
j=i sinh(i(4""^^ - 4"~^) - ^) sinh(i(4"~^^ + 4"~^^) - ??) 



n 



sinh(i(4""^^ - 4""'^) - sinh(i(4"'^) + 4"-^^) - 27?) 
t^k sinh(i(4"-^) - 4"-^)) + 2ry) sinh(i(4"~^) + 4"-^)) + 2i^) 
g sinh(u^"~^^ - uf^ + 2r/) sinh(4"~^^ + uf^ + 2r/) ^^^^ 
j=i sinh(ti[." — 4"^ — 2r/) sinh(n^"' ^'^ + 4"^ ~ 2r/) 



^ "gi sinh(4"^ - uf-^^ + 2j?) sinh(4''^ + uf'^^ + 2j?) 
j=i sinh(4"'^ — '"j""^'' ~ 2??) sinh(4"'* + 4""^^ ~ 2??) 

sinh(u^"''* — uf^ — 4??) sinh(u^"'' + 4"'' — 4r/) 
j^k sinh(n^ ' -Uj + 4rj) smh^u^ + ''^^ 

for d^^: 

^ _ sinh(l(4"-^^ - 4"-^^) + r?) sinh(l(4"-^^ + 4"-^^) + r?) 

j=i sinh(i(4"~^^ - 4""^) - ^) sinh(i(4''~^^ + 4""^^) " ^) 
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' smh(i(4' 



(n-1) (n-l)^ 



2rj) smh[^[Uf, + u 



2l^) 



(n-1) „,(n-l)^ 



n 



sinh(i(?74: 



(n-1) 



+ 2?]) sinh( 
2,7) sinh(l(4"-^) 



(n-1) (n-1) 



+ 2r/) 



2r?) 



sinh(i(4' 



('^-1) 



-2ry) smh(l(4"-^)+4")) 



-2r?) 



(30) 



^^^^ sinh(l(4") - 4^^-^)) + 2r?) sinh(i(4") 



+ 4"-')) + 2r/) 



^1 sinh(l(4")-4'^-^))-2r7) sinh(i(4") 



+ 4"-'^)-2r/) 



™" smh(i(u^"^ - 4"^) - 477) smh(i(4"^ + M^"^^ 



n • 1 / I/ (» 

jVfc sinh(2(M;fc 



4?7) 



) + 4r/) sinh(i(n^"^ + 4"^) + 4r?) 



(31) 



smh(i(4"' 


-2) 


(n- 
-4 


"■^^) +?]) sinh(i(4" 


-2) , (n- 
+ 4 


■3)) 


+ r]) 




-2) 


(n- 


sinh(i(4"~ 


■2) , (n- 
+ 4 


■3)) 


- v) 


smh(i(4"" 


-2) 


(n- 

-^• 


-'V2r/) sinh(i(4" 


-2) , (n 

+ n^- 


-2) 


)-2r/) 


■ 1 / 1 / 
sinli(2(4 


-2) 


(n- 


~'))+2r?) sinh(i(4" 


,-2) , (n 
+ 4 


-2) 


) + 2r/) 


sinh(2(Ufc 


-1) 


(n- 


-^))+r?) sinh(i(4"- 


■1) , (n- 
+ Uj 




+ r/) 




-1) 


(n- 
-4 


"^^) -77) smh(^(4"" 


■1) . (n- 
+ 4 







^ sinh(l(4-^) - uf) + r?) sinh(l(4"-^) + uf) + y) 
3=1 sinh(i(4""^^ - 4"^) - v) sinh(i(4""^^ + 4"^) - v) 



■ 1 / 1 / 
smh(2« 


-1) 


(n- 
-4 


"^^) +77) sinh(i(4" 


1) , (n-2)^ 




smh(i(4"" 


-1) 


(n^ 
-4 


^^^) -7]) sinh(i(4"~ 


1) , (ra-2)^ 


- v) 


smh(2« 


-1) 


(n- 
-4 


-^V277) smh(l(4" 


-1) , (n-1) 
+ 4 


)-2rj) 


■ 1 / 1 / 
sinh(2(^fe 


-1) 


(n- 
-4 


-^)) + 2r/) sinh(i(4" 


-1) , (n-1) 
+ 4 


) + 2r,) 



(33) 



-Q- »11i'Hi("A: - "j ) + ^"llHl ("/^' ' + "j + 

j=i sinh(i(4"^ - 4"~^^) - ^) smh(^(4"^ + 4""^^) - ^) 



9 



smh(i| 



n ,,1 



u 



27?) smh(i(4' 



(n) 



+ «S"^) 



27?) 



+ 27?) smh(i(4"^+^"0 + 27?) 



(34) 



One can therefore determine the unknown functions zi{u) 



sinh(u) sinh(u — 2k7?) cosh(u — u;7?) 
sinh(n — 2/7?) sinh(M — 2{l + 1)7?) cosh(u — kt?) 

/ = !,•• -,71-1 for 45_i 



sinh(n) sinh(n — 2^7/) sinh(ti — liJT?) 
sinh(ti — 2/77) sinh(u — 2{l + 1)7?) sinh(7i — k7?) 
/ = !,■• -,77-1 for 41), C«, Z)« 



(36) 



w{u) 



sinh(ti) sinh(ti — 2k7?) 
sinh(ti — 2777?) sinh(7t — 2(77 — 1)7?) 



for 



(37) 



We have obtained expressions for the transfer matrix eigenvalues, Eqs. (15) - (|2^), (|3q ) - (p7|) 



These expressions pass a number of checks which are similar to those performed in 
case. We are therefore confident that these are the correct eigenvalues. 



for the A. 



(2) 
2n 



3 

We conclude by listing some unsolved problems: 

The cases which remain to be treated are Dn \ D^l\ and (with the exception of ^2^'*) all the 
exceptional affine algebras. For these cases, the R and/or K matrices are not yet known. 

As noted in the introduction, the analytical Bethe Ansatz method has been further generalized 
III, ^ to certain open spin chains which are not quantum-algebra-invariant; namely, spin chains for 
which the K matrix is diagonal but is not necessarily equal to the identity matrix. (For a spin 
chain with a non-diagonal K matrix, the analytical Bethe Ansatz method presumably does not 
work, since an eigenstate (e.g., the pseudovacuum state) of the transfer matrix is not available.) It 
would be interesting to find new diagonal K matrices, and to diagonalize the corresponding transfer 
matrices. 

Other open problems include formulating the algebraic Bethe Ansatz for open spin chains (this 
is known only for the cases and B^^^); studying further examples of models with spins in 
higher-dimensional representations; and investigating "graded" models associated with solutions of 



10 



the graded Yang-Baxter equations. Perhaps the most interesting outstanding problem is to use the 
Bethe Ansatz results to investigate boundary phenomena in the thermodynamic (A^ oo) limit. 



We are grateful to M. Jimbo and A. Kuniba for valuable correspondence. This work was 
supported in part by the National Science Foundation under Grant PHY-92 09978. 



A The R matrix 

The R matrices associated with the fundamental representation of ^2n-i' Bn \ Cn \ Dn^ were 
found by Bazhanov |15| and Jimbo [p!^]. ^ We follow the latter reference; however, we use the 
variables u and ry instead of x and k, respectively, which are related as follows: 

X = , A; = e^" . (38) 

The R matrix is given by ^ 

R{u) = C{u) Eaa «> Eaa + b{u) ^ E^a ® ^/3/3 (39) 

+ (e(n) ^ +e(n) ^ )Eai3 ® Ej3a + Yaa(}{u)Eai3 ® E^' j3' 



with 



c(m) = 2sinh(f-2r/) 

b{u) = 2sinh(f) 

e{u) = —2e~2 sinh(2r7) 

e{u) = e"e(n) 



cosh(| 
sinh(| 



- Kl]) 
KT]) 



for A^S-i 

for B^r^\ C«, 



^(1) 



(40) 



^The A^^_-^ R matrix given in Jis] , is Uq{D„) invariant. We consider here a different ^2n-i ^ matrix, which 
instead is Uq{C„) invariant. We obtain |18| the latter R matrix from the R matrix by replacing (in the notation 
of the first paper in ^]) ^ — fc^"^^ by ^ = — fc^"; i.e., by changing ^ — » ~^k"^. The R matrix obtained in this 
way presumably coincides with the (7q(C„)-invariant A^2-i R matrix of Kuniba ^]. We remark that the A'>^^_'^ R 
matrix given in can be obtained from the one for o'il^ by changing ^ ^£,k^ ■ Similarly, the yljn ^ matrix, 
which is Uq[Bn) invariant, can be obtained from the B'^^ R matrix by changing ^ —^k^. 

^This expression for the R matrix differs from the one given in Ref. jig] by the overall factor 2e^'^^'^'^^^^ . 
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f o • / cosh(| - (AC - 2)r?) for 4n-i 

2sinh(^) X < V2 ^ fn L 

\ sinh(| - (At - 2)r?) for i^i^^ , , 



a = a 7^ a' 



6('u) — 2 sinh(2r7) sinh(2n — 1)?7 for 5, 



(1) 



a = I3,a = a! 



2sinh2rje^2 x < ' ^ «p vi '/ i 

I €„e^e(±«+2("-/3))'? sinh | - sinh(| - nri) for 5^^^ , C),> , D^n' 



where 



K = < 



2n for 

2n-l for B^^ 

2n + 2 for ci^' 

t 2n-2 for D^n^ 



(42) 



and 



a 



a + \ n + 1 < a < 2n for A^^l-i^ ci^^ 



(43) 



a = < 



a + i 1 < a < ^ 
a « = %i 



fori?«,Z^« 



(44) 



a, /3 



a 



iV + l-a 



(45) 



1 for 1 < a < n 
—1 for n + 1 < a <2n 



for 4^2) d'^^ 

lOr ^2n-l' '-^^ 



for B^^^ D^^ 



(46) 



where N = 2n for ^2n-ii ' ^'^'^ iV = 2n + 1 for Bn ; and the E^i} are elementary matrices. 
Evidently, the R matrix acts on the tensor product space ® . 

In addition to obeying the Yang-Baxter equation, this R matrix satisfies the following important 
properties: 
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PT symmetry 

Vi2 Ri2{u) Pi2 = R2i{u) = Ruiuf''' ; (47) 

unitarity 

Ri2{u) R2i{-u) = C{u), (48) 

where C(u) is given by 



^ ) -4sinh(|-2r/)cosh(|-K7?)sinh(f + 2r?)cosh(| + Kr/) for 4n-i ug) 
\ 4sinh(f - 2r?) sinh(f - nr]) sinh(f + 2r?) sinh(f + ktj) for Bi^\ dn \ ^ 

crossing symmetry 



Ruiu) = Vi Ri2{-u - pf^ Vi = Vi' Ri2{-u - pf^ , (50) 

where p = —iir — inrj for A2n_i and p = —2Kr] for Bn \ Cn \ Dn^; and V'^ = 1; 
regularity 

R{0) = -ao)'2V , (51) 

where V is the permutation operator 

r = Y,Ea(3^E/Sa; (52) 

R{u) , R{v)\ = 0, R = Vn- (53) 

ii:(n + 2m) = R{u) . (54) 

The crossing matrix is given by 

V = E^^S^^, - J2 e^"~""'^^^aa' • (55) 

Correspondingly, M = F is given by the N x N diagonal matrix 

M = c/ia5(e2(2")^,e2(2n-2),^... ^^4,^^-4,^... ^g-2(2n-2),^g-2(2n)^^ 

M = diag (e2(2«-i)'? , e^^^^-s)'? ,e^\l, e'^'' , • • • , e-2(2«-3)r; ^ g-2{2„-i)^^ ^(i) ^^g^ 

M = diag (e2(2'^-2)'? , e2(2n-4)'7 ,...,1,1,..., e-^(2n-i)r, ^ g-2(2n-2)„ j fo^ _ 
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